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METHOD AND SYSTEM FOR DETECTING SUBSTANCES, SUCH AS SPECIAL 

NUCLEAR MATERIALS 



FIELD OF THE INVENTION 
The present invention relates generally to methods and systems for detecting small 
amounts of special nuclear materials (SNM) hidden in cargo and vehicles. 



BACKGROUND OF THE INVENTION 

In recent years, ever-higher priority is being accorded to developing state-of-the- 
art security-screening techniques and devices that will reliably and rapidly detect small 
quantities (tens of grams) of SNM concealed within cargo and vehicles. 

The following are examples of the prior art: 

a) Passive radiation-monitoring systems: 

Such "Portal Monitors" were developed for the purpose of interdicting theft or 
diversion of SNM from high-security sites. The detection of SNM is based on measuring 
an increase in radiation intensity above the ambient background. Among others, LANL 
(Los Alamos National Laboratories) and LLNL (Lawrence Livermore National 
Laboratories) have been active in developing such systems, which are usually packaged in 
three distinct categories: 1) Small, Hand-Held Monitors 2) Automatic Pedestrian Monitors 
3) Automatic Vehicle Monitors. 

Vehicle monitoring systems for detecting SNM are now available from several 
commercial companies. Among these, Canberra's JPM-12A Vehicle Monitor was 
developed jointly with LANL. It consists of two large plastic scintillators positioned on 
either side of the portal, that measure the gamma-rays emitted from SNM. The instrument 
can detect 10 gr of Pu and 1000 gr of HEU in an unshielded vehicle moving at a speed of 
2 m/s. The disadvantage of this method is that gamma detection can be thwarted with 
relative ease by enveloping the SNM with Pb-sheet, since the minimum detected 
quantities increase rapidly by an order of magnitude when shielded by several mm of it. 
In the Pu case, the problem can be countered by using neutron detectors in the portal. 
However, due to low neutron emission rates, this approach requires large-area neutron 
detectors. LANL has developed such a system, based on two large-area 3 He detectors 
between which vehicles pass. An equivalent system is currently marketed by TSA- 
Systems-Ltd., Model NVM-245. NUCSAFE's Vehicle Monitor claims detection of 5 gr 
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of gamma-shielded weapons-grade Pu within 10 s at 1 m distance. Also available are 
gamma-neutron Vehicle Portal Monitors manufactured by LAURUS-Systems-Inc. Model 
No. VM-250GN and by POLIMASTER, Model Series PM5000. For the latter, which 
incorporate sizeable plastic-scintillator detectors as well as large-area 3 He proportional 
counters, detection of 4.3 gr 239 Pu and 300gr 23 *XJ at a scan speed of 10 km/hr is claimed. 
However, in analogy to y-rays, small Pu quantities might evade detection if some neutron 
shielding is introduced. 

b) Single-energy radiography systems 

This method yields information about the contents of the screened object. A 
number of systems have been developed in recent years by commercial companies. 
Nearly all of them use high energy bremsstrahlung radiation produced by high-power 
linacs, although SAIC have built a line of products (VACIS) around a radioactive 60 Co 
source. In several systems, two views at 90° to each other are generated, to obtain more 
information. Commercial companies such as Smiths-Heimann and Aracor produce linac- 
based vehicle and marine container inspection systems. They operate at electron energies 
of ~10 MeV and employ Cadmium-Tungstate radiation detectors. The spatial resolution 
permits the detection of 1 mm copper wire. The scan time is —3 minutes for a 20 m-long 
vehicle, but the visual inspection time of the resulting image typically takes 10-15 
minutes. These systems penetrate -30 cm of steel and the radiation dose to the object is 
about 150-250 p.Gy. Their obvious drawback is that they do not automatically identify 
SNM and in general, the performance of such systems relies heavily on operator skill and 
judgment. 

c) Dual-energy radiography (DER) systems 

DER is a well-established technique that has found numerous applications in 
medical imaging (in particular, in-vivo bone mineral densitometry), environmental 
studies, material assaying, NDT, NDE, as well as security inspection scenarios. They are 
based on comparing the transmission attenuation at two energies, a sensitive measure for 
the atomic number Z of the absorbers in the line-of-sight from radiation source to 
detector. 

DER scans are usually performed at energies where the photoelectric component 
of the transmission attenuation dominates. One particular variant exploits characteristic 
discontinuities in attenuation when the incident photon energy is varied around inner- 
atomic-shell (K,L,M,...,) binding energies of a particular element. Typically, the latter 
vary with Z . Locating such an energy-discontinuity (by comparing the photon flux 
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transmission above and below it) is a sensitive indication of the presence of the element in 
question in the field-of-view. However, this variant of the method is unsuitable for 
inspecting massive cargo items such as aviation or marine containers, since the .relevant 
K-binding energies (the highest-energy discontinuities) for Z ~ 92 are at -110 keV, an 
energy too low to penetrate the inspected items. 

In summary, the inadequacy of existing screening methods to effectively interdict 
the illicit-transport of SNM at the crucial outgoing and incoming control points 
underscores the need for novel inspection systems that will reliably detect the presence of 
small quantities of SNM in cargo items ranging in magnitude from small packages, 
through passenger baggage, palletized cargo and aviation containers, up to a full-size 
marine container, loaded land-vehicle or railroad freight-car. 

SUMMARY OF THE INVENTION 
The present work describes a novel and sensitive method for detecting small 
quantities of special nuclear materials (SNM) concealed in packages, passenger baggage, 
massive cargo and vehicles. It is based on dual-energy radiography (DER) of y-rays, 
whereby the pair-production-related increase in transmission attenuation with photon 
energy characteristic of high-Z materials is exploited to detect and distinguish them from 
the low-Z and medium-Z substances that make up the overwhelming majority of benign 
transported items. Moreover, it is possible to differentiate between threat materials and 
other, benign high-Z substances of lower density, such as Hg, Pb or Bi, provided at least 
two DER projections are taken. In order to construct a practicable system with good 
performance characteristics, judicious choices of y-ray energies and populating reactions 
employing low-energy ion beams at sub-rnA intensities are preferable. To this end, a 
comprehensive literature study of reaction-induced y-ray yields has been performed, 
providing best-mode radiation sources, although the invention is not limited to these 
radiation sources. These are the 12 C lines at 4.43 MeV & 15.09 MeV populated, in 
decreasing yield order, by the n B(d,n), 13 C( 3 He,a) & ,0 B( 3 He,p) reactions at energies 
Ebeam< 5 MeV . An alternative option is presented by the n B(p,y) reaction, for which E p 
-9 MeV may be needed to fully exploit the giant-dipole-resonance y-rays at 16-24 MeV. 
The yield of the latter reaction is lower than that of d-induced and 3 He-induced reactions, 
but it may give rise to superior image contrast sensitivity, the upper DER gamma-ray 
energy being higher. The best-suited detectors may be organic scintillators with pulse- 
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shape-discrimination (PSD) properties. As shown by a preliminary simulation, an 
operational DER inspection system built around such components may enable reliable 
detection of small (tens of gr) SNM quantities within short scanning times (typically, a 
few minutes), even in massive cargo items and loaded vehicles that exhibit high image 
clutter. Typical absorbed radiation doses to inspected items would be in the acceptable 
range of (1-10) juGy. 

There is thus provided in accordance with an embodiment of the present invention 
a method for detecting substances, the method including performing multi-view, multi- 
energy radiography by irradiating an object with a plurality of discrete high-energy 
gamma-rays at a plurality of different orientations, and detecting and mapping radiation 
passing through the object with at least one array of detectors, and indicating a presence 
of a high-Z substance (e.g., a special nuclear material (SNM)), by detecting a difference 
in a transmission attenuation characteristic of the high-Z substance as opposed to low-Z 
and medium-Z substances. 

The method may also include determining and localizing regions within the object 
containing the high-Z substance (or high density substance) with the multi-view, multi- 
energy radiography, in order to enhance specificity to SNM, as opposed to other high-Z 
materials 

There is also provided in accordance with an embodiment of the present invention 
a system for detecting substances, the system including a dual-energy radiography (DER) 
system including a gamma-ray radiation source, including a beam accelerator and a target 
to which the accelerator sends a beam thereby producing gamma rays, and a plurality of 
gamma ray detectors positioned to detect gamma ray beams that pass from the gamma-ray 
radiation source through an object to be inspected, wherein the DER system is adapted to 
indicate a presence of a high-Z substance, by detecting a difference in a transmission 
attenuation characteristic of the high-Z substance as opposed to low-Z and medium-Z 
substances. 

The DER system may make two measurements of transmission attenuation 
characteristics, one measurement performed at a photon energy close to the global 
absorption minimum for all atomic numbers Z (at ~4 MeV) and another at a higher 
photon energy. 
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BRIEF DESCRIPTION OF DRAWINGS 

In the drawings: 

Fig. 1 is a simplified graph of values of Aju as a function of Z for various photon 
energies, relative to 4.43 MeV, the abscissa being the atomic number Z, and the ordinate, 
Aju in cm 2 /g. 

Fig. 2 is a simplified graph of Ajj, • p , the difference in macroscopic attenuation 
coefficient calculated for Eyhigh =15.09 MeV and Eyiow = 4.43 MeV at the stable high-Z 
elements, beginning with the last rare earth (Hf), using nominal bulk densities, the 
abscissa being the atomic number Z, and the ordinate, Aju • p in 1/cm. 

Fig. 3 is a simplified graph of the systematics of ground-state Q-values for p- 
capture reactions on light nuclei, the abscissa being the target mass number A, and the 
ordinate, Q 0 in MeV. 

Fig. 4 is a simplified level scheme of l2 C showing the states of interest, along with 
various reactions that feed and de-excite them. The ll B(p,y) excitation curve is shown 
alongside the corresponding states, clearly showing the giant dipole resonances that decay 
to the ground-state and to the 1 st excited state. Also shown are the low-lying states of 12 B 
fed in n B(d,p) and their decay modes. 

Fig. 5 is a simplified graph of measured thick-target yields per projectile of 
reactions generating high-energy 12 C y-rays, the abscissa being E pr0 jectiie in MeV, and the 
ordinate, the yield (x 10" 9 ). The yellow (labeled Y) and magenta (labeled M) curves both 
pertain to the p-capture reaction, but the yields refer to gamma-rays leading to the 1 st 
excited state & ground state, respectively. The energies differ by 4.43 MeV and vary with 
incident proton energy as indicated near the curves themselves. The yield increase above 
~5 MeV is due to the giant dipole resonance. The blue (labeled B) and green (labeled G) 
curves refer to the 15.09 MeV yields from reactions with 3 He beams, of varying yield and 
spectral cleanliness, as indicated near the curves. The red curve (labeled R) represents the 
15.09 MeV yield from the (d,n) reaction on a U B target. Apart from the effects of 
neutrons, these gamma-spectra (see Figs. 6 and 7) are remarkably clean. For comparison, 
the black circle indicates the 9.17 MeV gamma-yield of p-capture on 13 C. 

Fig. 6 is a simplified graph of an in-beam Nal gamma-ray spectrum from 
bombardment of a U B target with a 2.75 MeV deuteron beam, the abscissa being the 
channel number and the ordinate, the relative intensity. 
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Fig. 7 is a simplified graph of an in-beam germanium gamma-ray spectrum from 
bombardment of a deuterium target with a 19.1 MeV n B-beam, the abscissa being Ey in 
MeV, and the ordinate, the number of counts. The equivalent d-beam energy is 3.47 MeV. 

Fig. 8 is a simplified schematic outline of a DER screening system, constructed 
and operative in accordance with an embodiment of the present invention, in the plane of 
the d-beam and central y-detector. 

DETAILED DESCRIPTION OF EMBODIMENTS 
Dual-Energy Radiography In The Present Scenario - Basic Physics 

In view of the fundamental limitations of existing screening equipment described 
in the previous section, recourse may be made to characteristic features of the attenuation 
process that come into play at higher photon energies (in the 1-30 MeV range), for which 
the radiation is much more penetrating. In this photon energy range, the pair-production 
attenuation mechanism comes into play, so that the total attenuation coefficient exhibits a 
rise above ~4 MeV, which is progressively more pronounced as the Z of the absorber 
increases. Thus, in order to distinguish between high-Z materials and lighter elements, 
one measurement may be performed at the global absorption minimum for all atomic 
numbers Z (at ~4 MeV) and another at higher energy. This is the basis for SNM detection 
underlying the present work. 

Specifically, the intensity-fraction transmitted through two absorbers of atomic 
number Z a , Z b at each of 'two different photon energies E y i & E y2 is given by the 
expression: 

T(E yi ) = e-^VPa-xa^b /P b.x b) ^ [=i2 eg (1) 

where p, jli & x are the density, specific attenuation coefficient and thickness in 
the line-of sight from source to detector, for absorbers a & b, respectively. The 
coefficients \x are characteristic and known functions of both Z and E y . 

It follows from eq. ( 1) that the Transmission Ratio at the two energies is given 

by: 

T(E yl ) / T(E y2 ) = e ~ (A ^ a * pa * Xa + A ^ b * pb ' Xb) eq. (2) 

Thus, high sensitivity for detecting absorber a (threat object) against a background 
of absorber b (benign object) will be obtained when the differences in attenuation 
coefficient Aju obey the following relations: 

I Ajj, a I = maximal and | Aja b I = minimal 
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or, alternatively, if: 

AjLi a , Aja b have opposite signs 

(clearly also, the effect of any benign absorber b on this ratio cannot be fully 
eliminated unless Afi b = 0). 

In order to illustrate the consequences of eq. (2) , the compiled values for the 
quantity: 

A|li = |Li(Eyhi g h) - p,(E y iow = 4.43 MeV) 
are shown in Fig. 1 as function of atomic number Z, for several photon energies 
Eyhigh. The choice of the lower energy E y iow = 4.43 MeV was dictated by the fact that it 

corresponds to the global minimum in attenuation coefficient p, as function of energy, and 
also to a copiously-produceable I2 C gamma-ray. 

Fig. 1 shows that the Ap, values are indeed small for medium-Z materials and even 
negative for low-Z materials. Moreover, the sensitivity to high-Z materials increases quite 

sharply with photon energy E y hi g h. Thus it is desirable, on contrast-sensitivity grounds, to 

employ the highest photon energy Eyhigh possible. 

In order to illustrate the specificity of DER to SNM, as distinct from other high-Z 
materials, Fig. 2 shows AMAC, the difference in macroscopic attenuation coefficient 
Ap- p (where p is the bulk density) of high-Z material traversed, calculated for E r hi g h = 
15.09 MeV and E y iow = 4.43 MeV. Such parametrization (see eq. (2)) is meaningful, since 
the volume distribution of AMAC for suspect items can be estimated with adequate 
precision (in voxels of ~1 cm 3 ) from a small number of DER views(at least two). 

Clearly, on the basis of the macroscopic attenuation, SNM (in Fig. 2, uranium and 
plutonium) can be readily distinguished from Rare-Earths and the heaviest stable 
elements (hafnium, thallium, lead and bismuth in Fig. 2) and to a lesser extent, from the 
transition metals (tantalum and mercury in Fig. 2). The only difficulty appears to be with 
the noble metals and their neighbours (tungsten, rhenium, osmium, iridium, platinum and 
gold in Fig. 2), but these are unlikely to be transported as undeclared cargo. 

Continuous-Spectrum And Discrete-Line Photon Sources For DER 

Hitherto, most DER applications have employed continuous-energy spectra, such 
as those produced in electron-beam-induced Bremsstrahlung (BS) sources. However, in 



8 

the present application, working with BS-sources is not preferred due to two of their 
fundamental characteristics, namely: 

I. the flux tends to be concentrated at low energies and tails off at high energies, 
thus reducing contrast 

IL the low duty factor of typical BS beams (~l%o) is unfavorable for spectroscopy 
(single event counting) 

Thus even though the present invention can be carried out with continuous- 
energy sources, nevertheless it is preferred to use discrete-energy nuclear-reaction- 
induced sources. Basic (although non-limiting) system requirements for DER are as 
follows: 

a. production of Eyiow and Eyhigh in the same nuclear reaction 

b. clean gamma spectra (in the sense of containing as few other lines as possible) 

c. minimal neutron production and activation of environment 

d. populating reaction accessible with low-energy, low-intensity, non-exotic beams 
and durable targets 

e. robust, reliable, compact and low-cost accelerator 

f. reliable, low-cost detectors that resolve Eyiow from Eyhigh and discriminate against 
interfering radiation 

These requirements on nuclear-reaction-induced sources, as well as the paramount 
contrast-sensitivity criterion, that favors working with the highest possible gamma-ray 

energy Eyhigh, serve as guidelines in selecting the most appropriate nuclear reactions for 

the present application. However, it is emphasized that the invention is not limited to 
these guidelines. 

Initially, the nucleus in which the optimal levels and transitions are to be found 
may be identified. This is readily accomplished by considering the proton separation 
energies in light nuclei, or equivalently, the ground-state Q-values (Q 0 ) for the (p,y) 
proton-capture reaction. The latter represent a measure of the intrinsic exo-thermicity of 
the nuclear process (over and above the kinetic energy imparted by the incoming proton) 
and thus of the degree of excitation achievable in the final nucleus. Their systematics in 
light nuclei, taken from a well-known compilation, are shown in Fig. 3. 
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The global variations and local fluctuations in this plot represent a superposition 
of nuclear-shell closure effects (at the "magic" numbers), pairing-force effects (odd-even 
staggering in neighboring nuclei) and most prominently, cc-clustering effects (the 
dominant repetition pattern of maxima in Qo is 4 mass units). 

The highest excitations (two open diamonds on the left in Fig. 3) are obtained for 
p-capture on 3 H and on 7 Li. However, both these reactions are not preferred for the 
present application, since the cross-section (and gamma-yield) of the former is extremely 
low and the latter leads to 8 Be, which breaks up into two a-particles without emitting 
gamma-rays. The next best target nucleus (by far) is 11 B (solid diamond just to the right of 
the above-mentioned open diamonds in Fig. 3). Thus the nucleus in which the highest 
excitations can be achieved is 12 C, which is the product of the n B(p ? y) reaction. 

Populating the final nucleus at high excitation is a necessary but not sufficient 
condition for generating high-energy gamma-rays with low-energy particle beams. A 
further requirement is that the highly-excited levels may decay predominantly (or at least 
appreciably) to the ground-state or to low-lying excited states in a single transition. In 
general, this is not the case, because, at high excitation, when a nucleus has the possibility 
of decaying by particle emission (via the so-called "strong" nuclear force), its branching 
ratio for emitting gamma-rays (decay via the "electromagnetic" force) is usually rather 
small. However, for certain "pathological" states, quantum-mechanical selection rules can 
drastically inhibit particle emission. Such states will indeed decay appreciably via gamma 
emission. Classic examples of this are isospin T = 1 states in self-conjugate (N=Z) light 
nuclei, which are forbidden to decay via oc-particle emission. Since T is known to be a 
good quantum number in light nuclei, these selection rules can inhibit particle emission 
by many orders of magnitude, and in certain cases they do. 

All these desirable properties are to be found in the I2 C nucleus (which has 
N=Z=6) 5 the level scheme of which is presented in Fig. 4. The high-energy gamma- 
transitions in 12 C, along with the 4.43 MeV transition, will constitute the first described 
embodiment of the present invention (see next section). Also shown in Fig. 4 are the low- 
lying levels of 12 B and their decay modes, which will figure in another embodiment, 
described further hereinbelow. 
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First Embodiment ~In-Beamu Reaction-Induced l2 C y-Ray Yields 

In the DER context, it is clear from Fig. 4 that the 4,43 MeV gamma-ray is a 

prime candidate for Eyiow . With respect to Eyhigh , the most attractive candidates are the 
discrete level at 15.09 MeV, accessible via the following reactions n B(p,y) , n B(d,n) , 
10 B( 3 He,p) and 13 C( 3 He,a) and the Giant Dipole Resonance (GDR) around 23 MeV 
(accessible via u B(p,y) at E p = 6 - 10 MeV). 

Nuclear Reaction Yields 

The expected thick-target gamma-ray yields (per projectile) of the 15.09 MeV line 
and the GDR lines in the above-mentioned reactions were obtained by energy-integration 
of the published experimental cross-sections and are shown in Fig. 5. 

Several comments on these reactions already appear in the caption to Fig. 5 and 
need not be repeated. Clearly, the choice of reaction defines an order-of-magnitude yield 

hierarchy: p-capture is the weakest (but the mean Eyhigh is higher than 15.09 MeV, 
thereby enhancing the contrast sensitivity) and 3 He-induced processes are intermediate. 
Deuteron-stripping is the best, although its feeding of the 15.09 MeV state has a 
threshold at E d -1.6 MeV. At E d -3 MeV, the yield of 15.09 MeV gamma-rays is ~3 
orders of magnitude higher than the 9.17 MeV yield in the Gamma-Resonance- 
Absorption (GRA) method underlying nitrogenous explosives detection. For those skilled 
in the art of security screening applications, this quantity probably best underscores the 
applicability of the present invention with the n B(d,n) reaction. At energies above E d ~3 
MeV, the thick-target (d,n) gamma yield will continue to rise approximately as E d K5 . 
This primarily reflects the deuteron range-dependence on energy, since the cross-sections 
are fairly flat. On the debit side, these reaction-induced gamma sources are characterized 
by relatively high associated fast-neutron yields, as illustrated by the following table: 



Table 1. Fast-neutron & 15.09 MeV y-ray Thick-Target Yields and their ratios at 
selected energies 



Projectile + 
Target 


Eprojectile 

[MeV] 


15.09 MeV y-ray Yield 


Fast-Neutron Yield 


Neutrons/ 
y(15.09 MeV) 


d + "B 


3.0 


4.5 • 10"* / d 


-2.3 • 10 s / d 


~5 


d + "B 


5.0 


6.0 • 10' 6 / d * 


<1.8 • 10" 5 / d * 


<3 * 
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3 He + 13 C 


5.0 


3.0-10" 6 / 3 He 


>210 s / 3 He ** 


>7 ** 


3 He + ,0 B 


5.3 


4.0-10" 7 / 3 He 


>4-10" 6 / 3 He** 


>10 ** 


p + n B 


3.0 


~1- 10' 8 / p (yo+Yi) 


0 (below (p,n) threshold) 


0 


p + n B 


3.5 


-3- 10 8 / p (yo+Y,) 


-1 • 10" 3 / p 


-300 


p + u B 


4.0 


-4- 10" 8 / p (Yo+Yi) 


-3 • 10 s / p 


-800 


p + n B 


6.0 


~1.6-10- 7 /p(Yo+Yi) 


-2.5 • 10' 4 / p 


-1500 


p + n B 


9.0 


7.0-10' 7 /p(Yo+Yi) 


-1.4 • 10" 3 / p 


-2000 


employing intermediately-thick target (~ 4 mg/cm"), which degrades d 


euteron beam 



energy down to 4 MeV 

** lower limits, since only ( 3 He,n) channel was taken into account, whereas ( 3 He,pn) was 
neglected (no data available) 

The neutron yields in Table 1 were calculated in analogous fashion to the gamma- 
ray yields of Fig. 5, namely, by energy-integration of the published experimental cross- 
sections. It is clear from the results that (d,n) has the edge over p-induced and 3 He- 
induced reactions not only with respect to gamma-yield, but also in terms of minimizing 
the neutron/gamma yield ratio. 

The unique exception to this trend is the p + 11 B system below the (p,n) threshold 
at E p =3.0 MeV, but the gamma yields per projectile are almost three orders of magnitude 
lower than with d + M B. With increasing proton energy, the gamma-yields increase quite 
sharply, but the neutron yields increase even faster, so that the neutron/gamma ratios 
become unacceptably high, as can be seen from the table. For 3 He-induced reactions, 
there is no bombarding energy at which the neutron/gamma yield ratio is minimized, 
since all gamma and neutron producing reactions are exothermic (no thresholds) and their 
cross-sections are essentially independent of energy. Finally, with the (d,n) reaction, there 
appears to be a possibility of reducing the neutron/gamma ratio without sacrificing 
gamma yield, by raising the beam energy to E d ~5 MeV and employing a target of a 
thickness that does not degrade the beam below E d ~4 MeV. This may be due to a marked 
reduction in the (d,n) cross-sections to the ground-state and first-excited-state of ,2 C in 
this energy range. Of course, such improvement might be at the expense of the 4.43 MeV 
gamma yield, which is undesirable. At higher deuteron energies the cross-section 
behavior is not known, so it is possible that even more favorable beam-target conditions 
might exist. However, this is unlikely, since at E d ~6 MeV the (d,2n) channel opens. 

The next topics to be considered are: a) whether the in-beam d + n B spectra are 
clean enough for the purposes of the present application and b) what measures can be 
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adopted to alleviate the adverse effects of the high neutron yield. These issues are 
discussed in the following paragraphs. 

Experimental Options for Ensuring Clean In-Beam Gamma Spectra 
In order to generate clean, background-free in-beam spectra, particularly with 
deuteron or 3 He beams, certain experimental measures commonly implemented in nuclear 
spectroscopy may have to be taken. In this context, it is instructive to examine two 
published, in-beam spectra. Fig. 6 displays a Nal gamma-ray spectrum with the n B(d,n) 
reaction at an energy of Ed - 2.75 MeV taken in the 1960's. Fig. 7 shows a spectrum 
taken in 2000 with "EUROBALL", a large Ge-spectrometer matrix surrounded by an 
array of Nal anti-coincidence counters. The latter employed the inverse reaction at a close 
center-of-mass energy. 

Some features of these spectra relevant to the present application are: 

a. the 15.09 MeV peak region is very clean in both cases 

b. the spectra are very similar for the two detectors 

c. the same characteristic continuum, extending up to -10 MeV equivalent y-ray 
energy, is seen 

d. the 4.43 MeV radiation is prominent 

e. in the Nal spectrum (Fig. 6), there are no other prominent peaks above 3.5 MeV 

The similarity of the Nal and Ge spectra, particularly with respect to the broad mid- 
spectrum bump, is a strong indication that this continuum is due to thermal-neutron 
capture in the detectors (which serve as "calorimeters" for the subsequent gamma 
cascades). The 3.09 MeV line in Fig. 6 was identified as being due to the 12 C(d 9 p) 13 C 
reaction, stemming from a carbon impurity in the boron target - indeed, as evident from 
Fig. 7, it is not seen in the inverse reaction. (Thus, care may be taken to prepare as pure a 
n B target as possible and to ensure that no carbon deposits accumulate on it during beam 
bombardment.) The source of the 5.5 MeV peak in the Ge spectrum (Fig. 7) is unknown, 
but this is not directly relevant to the problematics of the present case, since the inverse 
reaction will not be employed in this application. 

The above clearly shows that the primary sources of background under the 4.43 
MeV line are thermal-neutron-capture events and low-pulse-height response events to 
15.09 MeV gamma-rays. Measures that may be taken to alleviate the neutron-associated 
background problem are as follows: 
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1. using U B target thicknesses that do not degrade the beam below the 15.09 
MeV threshold (E d - 1-6 MeV) 

2. fast-neutron dumps up & downstream from target - u B(d,n) neutrons are 
forward & backward-peaked 

3. accordingly, locating the DER radiography system at angles around 90° to 
the deuteron beam 

4. surrounding the target by neutron-moderators and absorbers 

5. enveloping the detectors in a thermal-neutron-absorbing material, such as 
Li, B or Cd 

6. employing organic scintillators as gamma- ray detectors 

7. low-frequency pulsing of the deuteron beam (see the next embodiment) 

Measure #1 is explained in a paragraph further hereinbelow regarding targets - it 
should reduce neutron doses by -30%. 

Measures #2-5 are self-explanatory to the skilled artisan and will not be discussed 
further here. Each can be applied independently. 

Measure #6 is based on the favorable feature that low-Z materials, such as 
organic scintillators, do not exhibit the thermal-neutron-capture feature characteristic of 
high-Z detectors. In contrast, such detectors have the drawback of being efficient for 
detecting fast neutrons, primarily via proton recoil accompanying n-p scattering on 
hydrogen. Specifically, neutrons of energy around 10 MeV will produce pulses equivalent 
to -5 MeV gamma-rays, thereby generating spectral background under the 4.43 MeV 
region of interest. However, this problem may be solved by invoking the pulse-shape- 
discrimination (PSD) property of certain organic scintillators, that permits the distinction 
between fast-neutron and gamma-ray events in the detector via the difference in fall-time 
of the scintillation pulses. The PSD effect is most pronounced in liquids (where it gives 
rise to easily-detectable pulse-fall-time differences of 15-20 ns) but has also been 
observed, albeit to a lesser extent, in plastic scintillators. 

Once the fast and slow neutron-events have been suppressed, the principal spectral 
background under the 4.43 MeV line may be due to low-amplitude 15.09 MeV gamma- 
ray events. However, a response-optimized organic-scintillator may exhibit a flat 
background in the region of interest, that can be reliably subtracted from the 4.43 MeV 
peak. In this context, one should also consider the use of deuterated scintillators, in which 
the 2.2 MeV y-ray from slow-neutron-capture by hydrogen is eliminated. It would permit 
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utilizing a higher fraction of the response curve to 4.43 MeV gamma-rays, thereby 
increasing the detector efficiency. 

Employing organic scintillators might thus be advantageous (and also 
considerably cheaper), but may incur some loss in detection efficiency. 

It is noted that different kinds of scintillators may be employed in the invention. 
For example, as mentioned before, the scintillators may be organic scintillators (e.g., 
organic liquid scintillators and/or organic plastic scintillators, with or without pulse- 
shape-discrimination properties. Alternatively, the detectors used may include inorganic 
scintillator spectrometers, such as but not limited to, Nal, BGO or BaF2. As another 
alternative, the detectors used may include solid-state radiation spectrometers, such as but 
not limited to, Ge, CdTe or CdZnTe (CZT). When employing organic scintillators without 
the pulse-shape discrimination feature, or alternatively, when using inorganic scintillators 
or solid-state radiation spectrometers, invoking one or more of measures #1-5 may in 
itself provide adequate discrimination against spectral contamination by fast or slow 
neutrons. This depends on the detailed experimental conditions prevalent in a specific 
application of the invention. 

Implications Of y-Source For System Operating-Mode & Performance 

For the purpose of the discussion, a 2-dimensional schematic of the DER system 
embodiment is shown in Fig. 8. Its conventional aspects will not be elaborated on here, 
the emphasis being on features specific to this application. Fan beams and associated 
detector arrays are perpendicular to the plane of the drawing. 

Fig. 8 explicitly cites the n B(d,n) reaction, but would essentially be unchanged for 
p-capture or 3 He reactions. Several features will now be discussed (accelerator and target 
issues are discussed further hereinbelow). 

Multi-view DER 

One of the features illustrated in Fig. 8 is that, using a number of collimator and 
detector arrays, a set of radiographic projections (views) of the screened object may be 
obtained simultaneously. This translates into a gain in system throughput, as well as 
facilitating the density reconstruction procedure, which permits discrimination against 
other, innocuous high-Z substances (see also hereinbelow). 

Stringent collimation of the y-beams to a fan geometry (perpendicular to the X-Y 
plane of Fig. 8) may be used to minimize the adverse effects of scattered radiation, a 
major problem with high-energy gamma-rays, especially traversing large, massive 
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objects. Such scattering will primarily tend to reduce contrast sensitivity, but it will also 
have the effect of blurring the image features (edges). 
The PER object-screening sequence : 

With the interrogating radiation produced by a nuclear reaction (as opposed to a 
Bremsstrahlung source), the low amounts of radiation available may be judiciously 
utilized, if the DER system is to achieve high throughput, high detection probability and 
low false-alarm rates. 

The transmission sampling of the two y-ray energies may be effected in identical 
irradiation geometry through any part of the interrogated object. Accordingly, the object- 
screening sequence may comprise one or more (in most cases, not more than two) of the 
following stages: 

I. a rapid pre-scan, to classify the overall y-ray-attenuation general range 

II. a full scan with one detector array, to locate regions suspected of containing 
high-Z materials 

III. multi-view scans (Fig. 8), to establish presence of high-Z (along the lines 
previously described) 

IV. a localized scan with ultra-clean spectra (as mentioned above), to confirm 
presence of high-Z 

V. a repeat of one or more of the above stages in borderline cases, to eliminate 
persistent false-alarms 

At stages I - III, the object will be scanned in continuous translational motion, 
possibly at variable speeds. Stage IV will only be implemented on the regions that proved 
problematic in the continuous scans. 

In the DER system of the present invention, it may be feasible to subject all 
decision-making to computer control, with respect to: the specific scanning sequence to 
be performed, the mechanical motions required to execute it, as well as the data 
acquisition modes and analysis routines implemented for appraisal of the container status 
("as yet undetermined", "clean", "suspect" or "bingo") throughout its interrogation. 

Proiected performance characteristics : 

By way of illustration, the projected performance characteristics are estimated for 
the u B(d,n) reaction. 
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In-beam Gamma-ray yields from 

As can be read off Fig. 5, the 15.09 MeV y-yield / deuteron at E d = 3 MeV is 
|4.4 ■ 10*^ . In this context, it is noteworthy that the anisotropy of the in-beam 

15.09 MeV gamma-rays has been measured and found to exhibit very slight peaking 
(-10%) at 90° to the beam. 

Thus, for a 200 \xA, 3 MeV d-beam on a thick n B target, the 15.09 MeV yield is 

[5^10* 7s[ . 

The 4.43 MeV yields are slightly more difficult to estimate, since the 12 C first 
excited state is populated directly in the (d,n) reaction with a typical 25 mb cross-section 
(similar to that for the 15.09 MeV state), but is also fed by gamma cascades from higher 
states. Judging by the spectra in Figs. 6 and 7 and the limited cross-section data available 
for the higher states, the 4.43 MeV yield can be conservatively estimated as being the 
same as the 15.09 MeV yield (in reality, it may be a factor of 2-3 higher). 

Detector counting rates at full transmission {no absorber) 

The counting rates are estimated assuming organic scintillator detector dimensions 
of 1.5 x 1.5x40 cm 3 . 

Detector solid angle (1.5 x 1 .5 cm 2 cross-section) at 400 cm from target: 1,1 • 10 

Total intrinsic efficiency (40 cm length): 
at 4.43 MeV 

assumed useful response fraction: 
at 4.43 MeV 

multiplying the latter two entries gives: 
Useful intrinsic efficiency : 
at 4.43 MeV 



0.50 at 15.09 MeV and 



0.69 



!y$ at 15.09 MeV and h 0-2| 



15% at 15.09 MeV and 



14% 



Finally, taking yield x detector-solid-angle x useful-intrinsic-efficiency 



gives 



Full-transmission counting rates of : [880 /s| at 15.09 MeV and |810 /s| 

at 4.43 MeV 

In order to test the performance of the system in this counting-rate regime, 
detection and identification algorithms along the lines developed for GRA have been 
applied to synthetic images of an LD-3 aviation container and the expected probabilities 
of detection and false-alarm were determined. The simulations were performed for two 
different cargo contents: 
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Cargo Content # I 

LD-3 container, uniformly filled with low-Z (organic or aqueous cargo) 
equivalent to traversing 1 15 cm of water at every point in the scan (broadly representative 
of bulk agricultural produce, oil barrels, chemicals, etc.). 

Cargo Content # II 

An iron slab of dimensions 40x15x15 cm 3 , weighing -70 kg (representative of a 
heavy mechanical part). 

Into each cargo type, a cube of SNM (1.5x1.5x1.5 cm 3 in dimensions) weighing 
65 gr was introduced. 

For both these cases, 500 simulated DER scans of the container (each comprising 
two views at 90° orientation to each other) were generated using the counting- 
statistics/pixel corresponding to the calculated 4.43 MeV and 15.09 MeV transmission 
attenuations through the absorbers in question. 

The preliminary analysis performed here follows a procedure along the lines of 
the one developed for explosives detection via GRA. The estimated DER performance 
characteristics presented here are based solely on the detection and localization of the 
suspicious object. Table 2 below presents the results of this simulation. 



Table 2: Simulation results 



Cargo 
Content 


Assumed 
Counts/pixel 


Detection Probability 


False-Alarm 
Rate 


Scan- 
Time 


#1 


N 15 =120, N 4 =30 


>95% 


<1% 


~J min 


#11 


N 15 =150, N 4 =150 


>99% 


<1% 


—4 min 



where 7V/ 5 and N 4 are the assumed number of counts/pixel at the two energies 
detected after attenuation. The Scan-Time depends on the values assumed for A^/ 5 and N 4 , 
on the number of slices (-130) that constitute a full scan of the entire container length and 
on the unattenuated counting rates quoted above. 

It may be possible to enhance specificity to SNM by invoking criteria based on 
the object density, which is much higher for SNM than for high-Z materials such as Pb or 
Bi (see Fig. 2 and description above). To achieve such discrimination over voxels of 
typically 1cm 3 , a small number of radiographic views (2-6 or so) of the inspected object 
will suffice. Following Fig. 8, these could be taken in a single scan with multiple fan 
beams and detector arrays. It is thus possible to detect even smaller quantities of SNM 
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than the 65 gr cited above. However, even without this feature, the performance 
characteristics derived from the simulation are already very promising. 

Second Embodiment - Out-Of-Beam !2 C(4.43 MeV) Activation Yields 

The second embodiment describes in detail Measure #7 mentioned hereinabove, 
which is an option for scanning limited container regions with particularly clean gamma 
spectra, making use of in-beam prompt 15.09 MeV and out-of-beam delayed 4.43 MeV y- 
rays. This option is based on the same n B + d projectile-target system of the first 
embodiment, but via two different reaction channels. 

The following is an explanation of the principal feeding and decay modes 
underlying the spectral lines observed in Figs. 6 and 7. 

Apart from elastic scattering, the dominant nuclear reaction channels at incident 
energies E d < -5 MeV are: 

n B(d,n) 12 C and n B(d,p) 12 B for which T m = 20.2 ms 

The cross-sections for both (d,n) and (d,p) reaction channels are typically 20-30 
mb to each of the energetically-accessible states (including the 12 C 15.09 MeV state). 
They account for the intensities of the 0.95 & 1.67 MeV ( 12 B), 4.43 & 15.09 MeV ( I2 C) 
lines seen in the in-beam spectra of Figs. 6 and 7. One notable exception, however, is the 
cross-section for (d,p) to the 12 B ground-state. There is a major discrepancy regarding its 
value, as appears in prior art literature, ranging from ~25 mb in the E d = 1.0-2.6 MeV 
range, up to values 10-25 times higher (250 - 600 mb) at E d = 2.6 MeV. Another 
compilation, the EXFOR (CSISRS) database, quotes a cross-section of 176 mb at E d = 3 
MeV, decreasing slowly with increasing deuteron energy. Since all excited 12 B states 
decay to the ground-state within a fraction of a picosecond, the cumulative ,2 B production 
cross-section at each incident energy is the sum of the cross-sections to the individual 
excited states and the ground state. At E d =3 MeV, four such excited states are 
energetically accessible (see Fig. 4). Assuming 25 mb to each excited state, the 
cumulative l2 B production cross-section is thus, within the Aground-state cross-section 
discrepancy limits, in the range 100 - 700 mb. 

It is noted that the l2 B ground state decays to the 12 C(4.43 MeV) state with a 
branch intensity shown in Fig. 4 as 1.3% (the recently-adopted value is 1.23 (5)%). Thus, 
the u B(d,p) 12 B reaction can serve as a generator for 4.43 MeV activation gamma-rays, 
which may continue to be produced for a few tens of ms after the deuteron beam is 
switched off, until all 12 B nuclei have decayed. These delayed 4.43 MeV gamma-rays are 
the basis for Measure #7 above. 
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The value of the 12 B production cross-section may determine to a large extent 
whether the 4.43 MeV yields are sufficient for the present application. 

To illustrate this point, Table 3 shows the number of in-beam 15.09 MeV and out- 
of-beam 4.43 MeV activation gamma-rays produced in a thick n B target by a 200 pA 
beam of 3.0 MeV deuterons. They are respectively calculated for a 20 ms beam-burst, 
followed by 40 ms with no beam on target, assuming three different values for the 12 B 
production cross-section, namely: 100 mb (lower limit - column II), 250 mb (intermediate 
value - column III) & 700 mb (upper limit - column IV). 

Table 3: Number of gamma-rays produced in burst (parameter - 12 B 
production cross-section) 





cr P rod( 12 B) = 100mb 


a prO d( 12 B) = 250 mb 


c5 P rod( 12 B) =700mb 


Prompt 15.09 MeV: 
beam on (20ms) 


Delayed 4.43 MeV: 
beam off (40ms) 


Delayed 4.43 MeV: 
beam off (40ms) 


Delayed 4.43 MeV: 
beam off (40ms) 


1.1 • 10 8 


3.2 • 10 6 


8.0 • 10 6 


2.2 • 10 7 



Clearly, the out-of-beam spectra will be very clean, because there are no delayed 
gamma-rays at energies above 4.43 MeV. Thus, invoking measure #7 will furnish an 
option for measuring both relevant gamma-rays under background-free conditions - the 
15.09 MeV line in-beam and the 4.43 MeV line out of beam. The latter are obviously 
isotropic, being emitted from unoriented nuclei. 

However, as evident from Table 3, the out-of-beam yield of the 4.43 MeV 
gamma-rays is at least a factor of 5 (and maybe even a factor of up to -35, depending on 
the I2 B production cross-section) lower than that of the in-beam 15.09 MeV gamma-rays 
emitted during a single 20 ms beam burst. 



Counting-rate estimate: out-of-beam 4.43 MeV gamma-rays (in analogy to the 
section on in-beam gamma-ray yields from u B(d,nV ) 

Assume a prod ( 12 B) = 250 mb y No. of 4.43 MeV y-rays in 40 ms for 200 jliA of 3 



MeVd-beam = 8-10 e 



Detector solid angle (1.5 x 1.5 cm 2 cross-section, 40 cm length) at 400 cm 
from target = |1.1-10 6 j 
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Useful intrinsic efficiency : (assuming useful response fraction of 30%)= |21% 



we obtain : 1 -8 counts per (20+40) ms burst, or |30 /s| 

Thus, within the uncertainty limits on a pro d ( l2 B), we estimate: 
Full-transmission out-of-beam 4.43 MeV counting rates of : |l2 - 84/i 

Thus, in view of the low counting rates expected for out-of-beam 4.43 MeV y- 
rays, measure #7 is likely to be exercised primarily if recourse to a localized scan with 
ultra-clean spectra (as mentioned above), to confirm the presence of high-Z is indeed 
required to resolve the issue. In such a case, the localized scan might be performed in 
discrete stepping motion of the container, the number of 60 ms (20 ms beam-on + 40 ms 
beam-off) cycles, or total time spent at each scan region being determined by the 
attenuation regime and counting statistics required. Ultimately, the entire process may be 
automated. 

Accelerator And Target Options 

Subject to the nuclear reaction chosen, one of the following beams may be used: 

a. deuterons: 1 00 - 300 jjA of 3 - 5 MeV 

b. 3 He: 3 00 - 1 000 |iA of 3 - 5 MeV 

c. protons: 500 - 1500 pA of 9-10 MeV 

d. mixed d + & H 2 + : 2000 - 5000 \xA of 3.5 MeV 

Options a, b & c refer to reactions discussed in the section on Nuclear Reaction 
Yields, hereinabove, and will not be referred to further here. 

Option d might provide a solution for a particular application that calls for 
detection of high-Z materials and nitrogenous explosives in a single system, or even in 
the very same scan. In the latter case, one could bombard a thin 13 C layer deposited on the 
surface of a thick n B target, with a mixed mass = 2 beam of H 2 + (majority ion) and 
deuterons (minority ion). The H2 + molecular ions, when they impinge on the 13 C layer, 
will break up into two almost parallel, equal-energy protons at the 1.75 MeV resonance, 
as required for a GRA explosives detection system. Having traversed the I3 C layer, the 
protons will enter the n B target, where they will contribute to high-energy gamma-ray 
production via the p-capture reaction (see Fig. 5). The deuteron component will traverse 
the I3 C layer (contributing little in the way of reaction yield, the 13 C being thin) and 
interact with the n B target via the (d,n) reaction described above. 
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Alternately, should it be of interest to configure the high-Z and explosives 
detection systems around the same accelerator without performing the scans 
simultaneously, a simpler variant of the above is feasible: one could accelerate mass 2 to 
3.5 MeV in the machine, switching between hydrogen and deuterium in the ion source 
(and also changing the target) as required by the nature of the application at any given 
moment. 

In terms of accelerator technology, options a, b and c may be realized with 
virtually any type of machine: cyclotron, radio-frequency-quadrupole (RFQ) or 
electrostatic (single-stage Van-De-Graaff or two-stage Tandem Van-De-Graaff). Beam 
current requirements are not excessive for any of these types. In contrast, option d would 
probably favor an RFQ, as GRA requires considerably higher beam currents. 

Concerning targets, the requirements of the application are by no means excessive. 
High-quality and purity n B layers (self-supporting, if necessary) have been produced by 
vacuum-deposition for many years. Moreover, beam-heating problems should not be 
severe, since boron has an extremely high melting point. As mentioned in the description 
of the first embodiment, one of the options for reducing the neutron background is to limit 
the target thickness so that, when the deuteron beam energy is degraded to below the 
15.09 MeV yield threshold (at E d = 1.6 MeV), the beam emerges from the target and does 
not produce neutrons at lower bombarding energies. In practical terms, this means that the 
U B layers may be made ~2 mg/cm 2 thinner than the full deuteron range in boron (~5.5 
mg/cm 2 , for E d = 3 MeV ). Doing so will suppress -30% of the neutrons one would get if 
the target were thick enough to stop the deuterons (Measure #1 in the section concerning 
experimental options for ensuring clean in-beam gamma spectra). 

The preference for the n B(d,n) reaction (as mentioned in the description of the 
first embodiment) stems from its having the highest yield. However, an ultimate decision 
as to the optimal reaction may depend not only on its yield, but also on issues such as but 
not limited to, contrast-sensitivity, neutron doses to screened items and environment, 
accelerator complexity and cost, detector cost, etc. The accelerator considerations might 
eventually favor the p-capture and 3 He-induced reactions over (d,n), although higher 
beam currents will be needed to compensate for the lower reaction yields. 
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CLAIMS 

What is claimed is: 

1 . A method for detecting substances, the method comprising: 

performing multi-view, mult i -energy radiography by irradiating an object with a 
plurality of discrete high-energy gamma-rays at a plurality of different orientations, and 
detecting and mapping radiation passing through the object with at least one array of 
detectors; and 

indicating a presence of a high-Z substance by detecting a difference in a 
transmission attenuation characteristic of the high-Z substance as opposed to low-Z and 
medium-Z substances. 

2. The method according to claim 1, further comprising determining and localizing 
regions within said object containing the high-Z substance with the multi-view, multi- 
energy radiography. 

3. The method according to claim 1 or claim 2, further comprising determining and 
localizing regions within said object containing high-density substances with the multi- 
view, multi-energy radiography. 

4. The method according to any of claims 1-3, wherein the high-Z substance 
comprises a special nuclear material (SNM). 

5. A system for detecting substances, the system comprising: 

a dual-energy radiography (DER) system comprising a gamma-ray radiation 
source, including a beam accelerator and a target to which said accelerator sends a beam 
thereby producing gamma rays, and a plurality of gamma ray detectors positioned to 
detect gamma ray beams that pass from the gamma-ray radiation source through an object 
to be inspected, wherein the DER system is adapted to indicate a presence of a high-Z 
substance, by detecting a difference in a transmission attenuation characteristic of the 
high-Z substance as opposed to low-Z and medium-Z substances. 

6. The system according to claim 5, wherein said DER system is adapted to make 
two measurements of transmission attenuation characteristics, one measurement 
performed at the global absorption minimum for all atomic numbers Z (at ~4 MeV photon 
energy) and another at a higher photon energy. 

7. The system according to claim 6, wherein the global absorption minimum for all 
atomic numbers Z is approximately at 4 MeV photon energy. 

8. The system according to any of claims 5-7, wherein the transmission attenuation 
characteristic of the high-Z substance comprises a macroscopic attenuation coefficient 
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equal to the product of the specific attenuation coefficient and the bulk density of the 
high-Z substance. 

9. The system according to any of claims 5-7, wherein said gamma-ray radiation 
source comprises a discrete-energy nuclear-reaction-induced source. 

10. The system according to claim 9, wherein said nuclear-reaction-induced source 
has highly-excited levels that decay predominantly to a low-lying excited state in a single 
transition. 

11. The system according to claim 9 or claim 10, wherein said nuclear-reaction- 
induced source has highly-excited levels that decay predominantly to the ground-state in a 
single transition. 

12. The system according to any of claims 5-11, wherein said gamma-ray radiation 
source comprises at least one of the target and beam-projectile combinations n B+d ? 
13 C+ 3 He and 10 B+ 3 He , giving rise to nuclear reactions at beam energies Ebean^- 5 MeV. 

13. The system according to any of claims 5-12, wherein said gamma-ray radiation 
source comprises the n B(p,y) reaction at a proton beam energy of E p <~9 MeV. 

14. The system according to any of claims 5-13, wherein said gamma-ray radiation 
source comprises the n B(d,n) reaction, with a target thickness that does not degrade the 
beam below thel5.09 MeV threshold, wherein said threshold is at E d ~ 1.6 MeV. 

15. The system according to any of claims 5-14, wherein said gamma-ray radiation 
source also emits n B(d 5 n), and n B(d,n) neutrons that are forward and backward-peaked. 

16. The system according to any of claims 5-15, wherein said DER system is 
positioned at angles around 90° to the beam. 

17. The system according to any of claims 5-16, wherein said target is surrounded by 
at least one of neutron-moderators and absorbers. 

18. The system according to any of claims 5-17, wherein the detectors are at least 
partially enveloped in a thermal-neutron-absorbing material. 

19. The system according to any of claims 5-18, wherein said detectors comprise 
organic scintillators. 

20. The system according to any of claims 5-19, wherein said detectors comprise 
organic liquid scintillators. 

21. The system according to any of claims 5-19, wherein said detectors comprise 
organic plastic scintillators 

22. The system according to any of claims 5-21 , wherein said detectors comprise 
scintillators with pulse-shape-discrimination properties. 
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23. The system according to any of claims 5-18, wherein said detectors comprise 
inorganic scintillator spectrometers. 

24. The system according to any of claims 5-18, wherein said detectors comprise 
solid-state radiation spectrometers. 

25. The system according to any of claims 5-24, wherein said beam comprises a low- 
frequency pulsed deuteron beam. 

26. The system according to any of claims 5-25, wherein said DER system emits a 
beam of mixed d + and H 2 + ions at 3.5 MeV projectile energy, and wherein the DER 
system is adapted to indicate a presence of the high-Z substance and nitrogenous 
explosives as a function of the transmission attenuation characteristics of the high-Z 
substance and the nitrogenous explosives. 

27. The system according to claim 26, wherein the DER system is adapted to indicate 
the presence of the high-Z substance and the nitrogenous explosives in the same beam 
scan. 

28. The system according to claim 27, wherein the DER system is adapted to indicate 
the presence of the high-Z substance and the nitrogenous explosives by bombarding a thin 
l3 C layer deposited on the surface of a thick U B target, with a mixed mass comprising a 
beam of H 2 + and deuterons. 

29. The system according to any of claims 5-27, wherein the DER system is adapted 
to switch between ion sources and to switch between targets. 

30. The method according to any of claims 1-4, and substantially as shown and 
described hereinabove. 

31. The system according to any of claims 5-29, and substantially as shown and 
described hereinabove. 




For the Applicant, 




David Klein, Patent Attorney 
Dekel Patent Ltd. 
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